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We derive the higher dimensional generalization of Penrose-Tod equation describing past horizon 
in Robinson-Trautman spacetimes with a cosmological constant and pure radiation. Existence of 
its solutions in D > 4 dimensions is proved using tools for nonlinear elliptic partial differential 
equations. We show that this horizon is naturally a trapping and a dynamical horizon. The findings 
generalize results from D — 4. 
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I. INTRODUCTION 

Robinson-Trautman spacetimes represent a class of ex- 
panding nontwisting and nonshearing solutions [H-Q de- 
scribing generalized black holes. Various aspects of this 
family in four dimensions have been studied in the last 
two decades. In particular, the existence, asymptotic 
behaviour and global structure of vacuum Robinson- 
Trautman spacetimes of type II with spherical topol- 
ogy were investigated, most recently in the works of 
Chrusciel and Singleton In these rigorous stud- 

ies, which were based on the analysis of solutions to 
the nonlinear Robinson-Trautman equation for generic, 
arbitrarily strong smooth initial data, the spacetimes 
were shown to exist globally for all positive retarded 
times, and to converge asymptotically to a correspond- 
ing Schwarzschild metric. Interestingly, extension across 
the "Schwarzschild-like" future event horizon can only be 
made with a finite order of smoothness. Subsequently, 
these results were generalized in [7J, |8| to the Robinson- 
Trautman vacuum spacetimes which admit a nonvanish- 
ing cosmological constant A. These cosmological solu- 
tions settle down exponentially fast to a Schwarzschild- 
(anti-)dc Sitter solution at large times u. Finally, the 
Chrusciel-Singleton analysis was extended to Robinson- 
Trautman spacetimes including matter, namely pure ra- 
diation Q . It was demonstrated that these solutions with 
pure radiation and a cosmological constant exist for any 
smooth initial data, and that they approach the spheri- 
cally symmetric Vaidya-(anti-)de Sitter metric. 

In 10]. Robinson Trautman spacetimes (containing 
aligned pure radiation and a cosmological constant A) 
were generalized to any dimension. The evolution is 
governed by a simpler equation in higher dimensions, 
contrary to the four-dimensional case where fourth or- 
der parabolic type Robinson-Trautman equation oc- 
curs. Also, the possible algebraic types were deter- 
mined. But still several interesting features deserve at- 
tention, the presence of horizons being among them. Sim- 
ilarly to four dimensions, higher-dimensional Robinson- 
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Trautman family of solutions contains several important 
special cases, e.g. Schwarzschild-Kottler-Tangherlini 
black holes and generalizations of the Vaidya metric. The 
study of higher-dimensional spacetimes and their features 
help to comprehend which properties survive the gener- 
alisation and which are closely tied to four dimensions, 
thus deepening the understanding of General Relativity. 
Lately, considerable interest in higher dimensions comes 
from outside of the purely relativistic community. 

Our concern here is to locate the past (white hole) hori- 
zon. In general dynamical situations this might be rather 
nontrivial since the obvious candidate - event horizon - 
is a global characteristic and therefore the full spacetime 
evolution is necessary in order to localize it. Therefore, 
over the past years different quasilocal characterizations 
of black hole boundary were developed. The most im- 
portant ones being apparent horizon [llT |. trap ping hori- 
zon [l2| and isolated or dynamical horizon [l3l . Il4j |. The 
basic local condition in the above mentioned horizon def- 
initions is effectively the same: these horizons are sliced 
by marginally trapped hypersurfaces with vanishing ex- 
pansion of outgoing (ingoing) null congruence orthogonal 
to the surface. Quasilocal horizons are frequently used 
in numerical relativity for locating the black holes or in 
black hole thermodynamics. 

For the vacuum four dimensional Robinson-Trautman 
solutions without cosmological constant the location 
of the horizon together with its general existence and 
uniqueness has been studied by Tod [iEj . Later, Chow 
and Lun [l6| analyzed some other useful properties of 
this horizon and made numerical study of both the hori- 
zon equation and Robinson-Trautman equation. These 
results were recently extended to nonvanishing cosmolog- 
ical constant 1171. 



II. ROBINSON-TRAUTMAN SPACETIME IN D 
DIMENSIONS 

Robinson-Trautman spacetimes (containing aligned 
pure radiation and a cosmological constant A) in any 
dimension were obtained in [lO] using the geometric con- 
ditions of the original articles about the four-dimensional 
version of the spacetime [l|, 0|- Namely, they required 
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the existence of a twistfree, shearfree and expanding null 
geodesic congruence. They have arrived at the following 
metric valid in higher dimensions 

r 2 

s 2 = Hj *V - 2 ur - 2H u 2 (2.1) 

where 2H = (D _ 2 f (D _ 3) -2r(lnP) itl - (U^z^T) r2 ~ 

. The unimodular spatial (D — 2)-dimensional met- 
ric Jij(x) and the function P{x,u) must satisfy the field 
equation TZij = jjz^hij (with hij = P~ 2 jij being the 
rescaled metric) and n(u) is a "mass function" (we as- 
sume fi > 0). In D = 4 the field equation is always 
satisfied and 1Z (Ricci scalar of the metric h) generally 
depends on x 1 . However, in D > 4 the dependence on 
x % is ruled out (1Z = lZ(u)). But generally it still allows 
a huge variety of possible spatial metrics (e.g., for 
1Z > and 5 < D — 2 < 9 an infinite number of com- 
pact Einstein spaces were classified). The dynamics is 
also different in D > 4. While in four dimensions there 
is a fourth order Robinson- Trautman equation, the cor- 
responding evolution equation is much easier in higher 
dimensions 

N „ „N 167m 2 

(D - 1) p (lnP), u - = — - , (2.2) 
where function n describes the aligned pure radiation. 



III. PAST HORIZON 




FIG. 1: Schematic conformal diagram of Robinson- Trautman 
spacetime in D = 4 with A = and indicated past (trapping) 
horizon (PH) and event horizon (EH). 



The explicit parametrization of the past horizon hy- 
persurface is r = R(u, x % ) such that its intersection with 
each u — u\ slice is an outer marginally past trapped 
(D — 2)-surface. 

For the calculation of the expansion of an appropri- 
ate null congruence we will use a straight-forward gen- 
eralization of the tetrad formalism to arbitrary dimen- 
sion. Note that one can no longer use complex vector 
notation. Using two null covectors l a ,n a (with normal- 
ization l a n a = — 1) and D — 2 spatial covectors m a ^ 
(i = 1, ..,D — 2) we suppose the following decomposition 
of the metric 



In our case, we will be dealing only with the condition 
of vanishing expansion defining the marginally trapped 
hypersurfaces. Concretely, we will search for the past 
horizon similarly to previous studies in four dimensions 
and corresponding to the form of the metric containing 
retarded time. As will be clear later, we might call it trap- 
ping horizon or even dynamical horizon if it is spacelike, 
assuming appropriate higher-dimensional generalization 
of these notions (see [13|)- In four-dimensional case 
the parabolic character of Robinson-Trautman equation 
makes it generally impossible to extend the spacetime 
to past null infinity (the solutions of the Robinson- 
Trautman equation are generally diverging when ap- 
proaching u = —00) and it is impossible to define 
event horizon. In higher dimensions this is no longer 
truth (the evolution equation is different) but since one 
would like to investigate the horizon existence generi- 
cally, without prior specification of all necessary func- 
tions (e.g. dynamics of pure radiation) and geometry of 
(D — 2)-dimensional spatial hypersurfaces, the best ap- 
proach is still using the quasilocal horizons. In figure [TJ 
the schematic conformal picture of Robinson-Trautman 
spacetime (for D — 4 and without cosmological constant 
for simplicity) is presented together with the approxi- 
mate location of the horizons (initial data are given at 

U = Uq). 



g a b = -2 l( a n b ) + m a{i} m b{j} <P J (3.1) 

Null D-ad (D-bein) adapted to the trapped hypersurface 
(using the above mentioned parametrization) has the fol- 
lowing form: 

l a = (0,1,0, ..,0) 

n a = (l, [-H + ^h(VR,VR)} , VRj (3.2) 
m| i} - (0,Prh(yR,Wi),±Wi) 

where D — 2 vectors Wj diagonalize metric h, Vi? = 
{R' x ,.-,R ,x } and h(-,-) denotes scalar product 
w.r.t. h. Fortunately, in subsequent calculations we do 
not need the explicit form of the vectors w s , it is sufficient 
to know their orthogonality properties. 

By straight-forward computation one easily calculates 
the expansion associated with the congruence generated 
by l a to be 9; = ^f^- meaning that the outgoing null 
congruence is diverging. This is exactly what one as- 
sumes when dealing with the past trapped surface and 
is the additional condition in the definition of trapping 
horizon [12j. 
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IV. GENERALIZED PENROSE-TOD 
EQUATION 

Ingoing null congruence expansion can be calculated 
using the formula (sometimes a (D — 2) factor is used in 
the definition, but we are going to evaluate it to zero any- 
way) 9„ = n a -bP ab , where the tensor p ab = g ab + 2 l^n^ 
corresponds to the hypersurface projector. From 0„ = 
(equivalent to Penrose-Tod equation in four dimensions) 
we get the marginally trapped hypersurface condition 

n _ 2<g-3) ar2 _ (D _ 2){D _ 3) _V__ 2{D _ 3)A (i nR y 

- (D-4){D-3)h(VlnR,VlnR) = (4.1) 

It is a nonlinear second order partial differential equation, 
where both the laplacian and scalar product in the last 
term correspond to the Einstein metric h%j. Interesting 
property of this equation is that for D > 4 its nonlinear- 
ity is much worse since the term quadratic in derivatives 
appears. 

A. Results for D — A 

In four-dimensional case one can no longer use the ex- 
istence proof for equation (|4.1[) given by Tod il5| when 
the cosmological constant is present. However, one can 
use the version of sub and super-solution method adapted 
to Riemannian manifolds by Isenberg [l!| and valid for 
equations of the form A?/> = f(x,ip). For the proof 
of uniqueness one may use a straightforward modifica- 
tion of the original Tod's proof [l5| incorporating the 
cosmological constant. Using Newmann-Penrose equa- 
tions one can also determine the character of the horizon 
as a three-dimensional hypersurface. These results (for 
u = 2m = const.) are derived in |17| and summarized in 
the following table: 

TABLE I: D = 4 



RESULTS 


A = 


A < 


A > 


Existence 


Always 


Always 




Uniqueness 


Always 


Always 




Spacelike or null 


Always 


Always 





The restrictions for the positive cosmological constant 
can be easily understood by specializing to spherical sym- 
metry and A > (Schwarzschild de-Sitter) : 

• A < TT^- = rules out an over-extreme case. 

• R < ylfc = for the extreme case (9Am 2 

1 . reduces to R < 3m which may be interpreted 



as showing the uniqueness of the past black/white 
hole horizon (as opposed to the cosmological one). 

Both explanations are quite natural and not surprising. 

B. D > 4 : Existence of the solution 

The methods used in D = 4 are not applicable when 
the equation is of the form (after the substitution R = 
Ce~ u in (|4.1j) . assuming u > with a suitable constant 
C) 

Au = F(x,u,Vu) , (4.2) 

where F is quadratic in gradient. 

To prove existence of the solution to this quasilinear 
equation we will proceed by combining several steps (mo- 
tivated by [13] and using results from [2l| - |23l |). 

1. We will consider the differential operator 
Pu = —2(D — 3)Au + pu, with p > on a 
Riemannian manifold M. By using Maximum 
Principle we can prove that ker(P) =0 21]. 

2. The linear differential equation Pu = f with 
/ e C°' a (M) (Holder space over M) has unique 
solution u G C 2,a (M) (this standard result can be 
proven for example by Fredholm alternative [2l| 
and the previous step). 

3. To proceed with the nonlinear problem Pu = 
f(x, u, Vu), with / determined from (|4.1[) as (|| • \\h 
stands for the norm with respect to the positive 
definite metric hij) 

f = -pu + TZ-^§^AC 2 e- 2u ~ 

- { D-2)(D-3)nC Z - D e( D -^ - {D-A){D-S)\\Wuf h , 

we introduce the following truncature [2(], [22[ : 
f n - truncature of / by ±n. 

Then the map v 6 C^M) -> f n (x,v,Vv) 
is bounded. Using the previous step together 
with results on composition of Holder functions 
there exists a unique w e C 2,Q/3 (M) solving 
Pw = fn(x,v,\7v). 

4. The map v — ¥ w from previous step satisfies con- 
ditions of Schauder Fixed Point theorem, namely 
the a priori boundedness (see [2l| or (23|) =>■ for 
each n there is a fixed point u n G C 1,,3 (Af) (even 
u n e C 2 < a P(M)) solving Pu n = f n (x, u n , Vu n ) and 
moreover one can easily verify that ||u n ||x,oo < ^ 
(considering Pu n — f n (x, v, Vv) < n and compact- 
ness for evaluation at the maximum of u n ). 
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Using results of Boccardo, Murat & Puel [22||, in 
particular their Proposition 3.6, we can state the 
following corollary: 

Assuming that metric hij is smooth, function F 
can be estimated like \F\ < B(u)(l + | Vu| 2 ) (where 
B(u) is increasing function on R + ), and there ex- 
ist a sub- and a super-solution [24| u~ < u + , 
it* £ C l 'P(M) n L°°(M), then there is a L°°- 
bounded subsequence of the approximating so- 
lutions from the previous step satisfying u~ < u fl < 
u + a.e. 

Indeed, inspecting the above defined function / = 
pu — 2(D — 3)F one can verify that function B(u) 
might be found, namely there is no singular be- 
haviour at u = 0. Also, the domain we are dealing 
with is compact and therefore any dependence on 
x can be bounded for well behaved objects we use. 
For example, one may select the following bounding 
function 



B(u) = max W, 



(£)-2)^C 3 



D-A 



max 

xeM 



\hij(x)\ 



+ 



,(D-3)u 



where W 



AC 
D-l 



(_D-2)pC 3 



matrix norm of h was used. 



and the 



> F 



2, mm 



|Vu fi 



where the left hand side might be independently 
estimated. According to previous results gradient 
of u„ is bounded and from the last equation even 
independently. Therefore, the constant N does not 
depend on n. 

7. Then there is a C 1, ^-convergent subsequence Un — > 
u s , which proves the existence of the solution pro- 
vided the sub- and super-solutions are obtained. 
Moreover, using the second step with f(x, u s , Vu s ) 
we must have u s £ C 2 '"^(M). 

As is most common in the literature, we would be look- 
ing for constant sub- and super-solutions, first in the case 
TZ > (assuming u m i n > which can always be arranged 
by a suitable choice of C) : 

• A < 



75^3 ln 



C L 



C L 



(D-2)(_D-3)ai 



n 



(£>-2)(D-3),u 



(TZ 



2(D-3) 
D-l 



AC 2 



A > 



6. Thanks to elliptic estimate ||it fi ||c2, T < if (||u ft ||(70 + 
WUWcon) <K(\\u + \\ c o+N\\f\\ c o„)itis evenC 1 ^- 
bounded. To estimate f„ in the last inequality one 
can use its representation as f Tl = fg n , where func- 
tion 



9n 



1 



e(f-n) (l-f)-6(-/-n)(l + f) 



is responsible for the truncation and is the Heav- 
iside function. Using the results for composition 
(e.g. Holder index of composed map is a prod- 
uct of indices of components) and multiplication 
(e.g. index is a minimum of indices of compo- 
nents) of Holder continuous functions on bounded 
sets there has to be a new Holder coefficient 7 < a/3 
and a suitable constant N fulfilling the inequal- 
ity. Function / in the elliptic estimate is depen- 
dent on x not only explicitly but also via 
and Vun(x), which is reflected in the constant N. 
While Un is bounded independently on n by 
we need to bound the gradient in the same way. 
Using the fact that our function F has the form 
Fi(u) + F 2 (x)\\Vu\\j l (with strictly positive F 2 ) and 
integrating over the manifold (using the Stokes the- 
orem to eliminate the laplacian) we get 



/ Fxiun) = / F 2 ( x )\\Vu n \\l > 

IM JM 



These solutions satisfy all the conditions for any A < 0, 
but for positive cosmological constant one has to demand 



2TZ / TZ \ 

(U-l)(r>-2)(_D-3)Ai V2Ay 



> 1 



(4.3) 



so that u 2 > is valid for such a constant C which 
is maximizing the value of u\ '. Interestingly, this last 
condition reduces in the four-dimensional case (that was 
not explicitly studied here, but can be included trivially 
- note that then TZ is not a constant ona = const, and 
u ± has to be adjusted |17| ) to the condition from the 
table U for the existence of the solution when A > 0. 
One has to remember that for D = 4 scalar curvature TZ 
asymptotically (as u — > 00) approaches value 2. 

One may wonder whether the condition (|4.3[) is nec- 
essary or if it might be weakened by the choice of more 
suitable non-constant sub- and/or super-solution. Let us 
assume we have a positive solution R E C 2 ' S (M) of (|4.I[) 
with A > 0. Since it represents a function on a compact 
manifold it has to attain its maximum Rmax and mini- 
mum Rmin- At R m in the gradient term in (|4.1[) vanishes 
while A(\nR m i n ) > leading to the following inequality 



TZ 



2{D - 3) AR 2 - 

D—l L mm 



+ (D_2)(D-3)Mifcf' < 



(4.4) 
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The right-hand side of (|4.4|) has minimum at R mi7 ^ E — 

(fl-l)(^-2)(D-3)^ inequality gj} must alsQ hold 

for this value and after its substitution one arrives ex- 
actly at the condition (j4.3[) . Therefore, it represents not 
only sufficient but also necessary condition for the exis- 
tence of the horizon when positive cosmological constant 
is present. 

Since according to mathematical results any manifold 
(including compact ones) of dimension greater than or 
equal to 3 can be endowed with a complete Rieman- 
nian metric of constant negative scalar curvature [25l . [26[ 
one should also consider that TZ < for our (D — 2)- 
dimensional spatial hypersurface. One can propose the 
following constant sub- and super-solutions for 1Z < 
(assuming u mm > 0) 



• A < 



2C 

' (D-l)(D-2)^ 



A 



and select C > C m m that is defined by 

n _ 2(13-3) x n 2 



o x A.CZ Un -(D-2)(,D-a)nCjg 3) =0 



• A > : impossible to find constant u + . 

So for nonpositive scalar curvature 1Z we can prove the 
existence only for negative cosmological constant. 

In the four-dimensional case one can infer some useful 
results from the previous inequalities, mainly due to the 
fact that one can bound the scalar curvature from the 
asymptotic behaviour or using Gauss-Bonnet theorem. 
However, in higher-dimensional spacetime neither tool is 
available (the generalizations of Gauss-Bonnet theorem 
are very complicated and not immediately applicable) . 



C. D > 4 : Character of the horizon 

After establishing the existence of the past horizon T-L 
as a hypersurface foliated by marginally trapped surfaces 
one is naturally interested in whether it satisfies other 
conditions of recent quasilocal horizon definitions. We 
will consider trapping and dynamical horizons. 

The previous results tell us that 9; > and 0„ = 
holds on the past horizon. Since the Lie derivative CiQ n 
is in general nonvanishing on the horizon its closure is a 
trapping horizon 12j. Moreover, we can try to determine 
whether £i@ n < on the horizon which would mean that 
it is outer trapping horizon. After simple manipulations 
one arrives at the following formula 



-2- A 



(D-3)(D-2) 



AR 



D-5 



R D - 3 
(4.5) 



and we need to prove that 



D-l 



AR 



D-3 



(D-3)(D-2) 
2 



LiR- 



< 



D-5 1 



(4.6) 



Integrating the last equation over the D — 2 dimensional 
compact subspace spanned by coordinates x l (thus elim- 
inating the right-hand side in (|4.6p ) we get the necessary 
condition for the horizon being outer 



D-l 



AR 



D-3 



(D-3)(D-2) 



HR~ 2 < 



(4.7) 



which is satisfied for any A < and for positive cosmo- 
logical constant one shall demand 



R 



D 



-1 (J>-l)(P-2)(£>-3) M 
^ 4 A 



(4.8) 



Alternatively, one can consider (I4.6[) at the maximum 
of R, where AR D ~ 5 < 0. Notice, that for non-negative 
cosmological constant the right-hand side of (|4.6I) is 
strictly increasing function of R, so it has maximal value 
at maximum of R. If the horizon is everywhere non- 
degenerate (Ci& n 0) then it is really an outer trapping 
horizon. 

Next, we will consider a gradient of the horizon hyper- 
surface which in our parametrization reduces to 



N = r - R. u u - R jX J 



(4.9) 



We can use the sign of its norm to determine the causal 
character of the horizon. Using the null D-ad (I3.2|) the 
corresponding vector can be expressed in the following 
form 



N=UN a N a )\-n 



(4.10) 



and beside being normal to the horizon H it is also or- 
thogonal to its u = const. D — 2 dimensional sections 
R u . One can introduce a second vector orthogonal to 
these sections 



Z = \(N a N a )\ + rv , 



(4.11) 



which satisfies N a Z a = and therefore is tangent to the 
horizon %. Then (inspired by [Hj]) £z@n\u = holds 
identically which gives the following equation 



^(N a N a )de n + c n e n = o 



(4.12) 



If the outer trapping horizon condition is satisfied 
(£|6 n < 0) we need to determine the sign of the second 
term £„O n . One can consider the higher-dimensional 
generalization of Raychaudhuri equation (27| which in 
the case of nontwisting and nonshearing solution on the 
horizon where O n = simplifies to 



C n ®n = -Rabn a n b 



(4.13) 



Since we are considering only aligned null radiation (in 
the direction of 1 with radiation density $) and cosmolog- 
ical constant A the Ricci tensor in (|4.13p can be written 
as R a b = Ag a b + &l a lb and substituting into (|4.13[) we get 



C O 



(4.14) 
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Therefore, for nonnegative radiation density $ we con- 
clude that £„6 n < 0. Since both Lie derivatives in (I4.12[) 
are nonpositive it follows that N a N a < which means 
that the horizon is either null (for $ = 0) or spacelike (for 
$ > 0). In the latter case it presents an explicit example 
of dynamical horizon in higher-dimensional spacetime. 

V. CONCLUSION AND FINAL REMARKS 

We have derived the generalization of the Penrose- 
Tod equation to higher dimensional Robinson- Trautman 
spacetimes including cosmological constant and pure ra- 
diation. Using several mathematical tools we have proved 
the existence of its solution for any A < and for 
1Z > 0. The limitations arising for positive A (|4.3[) are 
shown to correspond to similar restriction arising in four- 
dimensional case that are naturally related to the more 
complicated horizon structure of relevant spacetimes (e.g. 
naked singularities) . Since the sign of scalar curvature of 
(D — 2)-dimcnsional spatial hypersurface does not restrict 
their topology as it does in D = 4 we have included the 
nonpositive case as well. 

Additionally, we have proved that one can consider 
this horizon as being a higher-dimensional generalization 
of trapping and dynamical horizon provided additional 
conditions are satisfied. 

The results show that in terms of the presence of the 
quasilocal horizons the higher-dimensional generalization 
shares the same qualitative behavior as the standard 
four-dimensional Robinson- Trautman spacetime. This 



provides support for considering the generalization given 
in (lpj to be natural not only mathematically but also 
physically. 

Several important issues were not investigated here, 
namely uniqueness of the horizon hypersurface and its 
possible topologies. The question of uniqueness is much 
harder to solve for D > 4 because of the nonlinearity 
in gradient. Due to our parametrization of the horizon 
the issue of its topology is connected with the topology 
of the underlying spatial geometry (given by hij) of the 
(D — 2)-dimcnsional manifold M. So the obvious start- 
ing point should be the classification of Einstein spaces 
of corresponding dimension. For D = 5 the S 1 x S 2 
(black ring) is ruled out since it cannot be endowed with 
Einstein metric [21| (the second homotopy group has to 
vanish ^(M) = 0) and the Poincare conjecture singles 
out three-sphere as the only simply connected case. In 
D = 6 topological obstructions arise for example due 
to generalized Gauss-Bonnet theorem relating the Eulcr 
characteristic x{M) and curvature of compact oriented 
four- manifold. It turns out that x(M) > and it is zero 
only in the flat case. This rules out S 1 x S 3 . In higher di- 
mensions the restrictions are much weaker (positive Ricci 
curvature implies finite first homotopy group). 

Acknowledgments 



This work was supported by grant GACR 202/09/0772 
and the Czech Ministry of Education project Center of 
Theoretical Astrophysics LC06014. 



[1] I. Robinson and A. Trautman, Phys. Rev. Lett. 4, 431 
(1960). 

[2] I. Robinson and A. Trautman, Proc. Roy. Soc. Lond. 

A265 , 463 (1962). 
[3] H. Stephani, D. Kramer, M.A.H. MacCallum, C. Hoense- 

laers and E. Herlt, Exact Solutions of the Einstein's Field 

Equations, 2nd edn (CUPress Cambridge, 2002). 
[4] P.T. Chrusciel, Commun. Math. Phys. 137, 289 (1991). 
[5] P.T. Chrusciel, Proc. Roy. Soc. Lond. A436, 299 (1992). 
[6] P.T. Chrusciel and D.B. Singleton, Commun. Math. 

Phys. 147, 137 (1992). 
[7] J. Bicak and J. Podolsky, Phys. Rev. D 52, 887 (1995). 
[8] J. Bicak and J. Podolsky, Phys. Rev. D 55, 1985 (1997). 
[9] J. Podolsky and O. Svftek, Phys. Rev. D 71, 124001 

(2005). 

[10] J. Podolsky and M. Ortaggio, Class. Quant. Grav. 23, 
5785 (2006). 

[11] S.W. Hawking and G.F.R. Ellis, The large scale structure 
of space-time (CUPress Cambridge, 1975). 

[12] S.A. Hayward, Phys. Rev. D 49, 6467 (1994). 

[13] A. Ashtekar , C. Beetle and S. Fairhurst, Class. Quant. 
Grav. 17, 253 (2000). 

[14] A. Ashtekar and B. Krishnan, Phys. Rev. D 68, 104030 



(2003). 

[15] K.P. Tod, Class. Quantum Grav. 6, 1159 (1989). 
[16] E.W.M. Chow and A.W.C. Lun, J. Austr. Math. Soc. B 
41, 217 (1999). 

[17] J. Podolsky and O. Svftek, Phys. Rev. D 80, 124042 
(2009). 

[18] L.M. Cao, JHEP 3, 112 (2011). 

[19] J. Isenberg, Class. Quantum Grav. 12, 2249 (1995). 

[20] T-H. Kuo, Nonlinear Anal. 63, e427 (2005). 

[21] A. Besse, Einstein Manifolds (Springer Berlin, 1987). 

[22] L. Boccardo, F. Murat and J. P. Puel, Annali della Scuola 
Normale Superiore di Pisa (Classe di Scienze) Ser. 4 vol. 
11 no. 2, 213 (1984). 

[23] D. Gilbarg and N.S. Trudinger, Elliptic Partial Differen- 
tial Equations of Second Order (Springer Berlin, 2001). 

[24] E.g. supersolution satisfies Pu + < f(x, u + , Vu + ). 

[25] T. Aubin, J. Funct. Anal. 32, 148 (1976). 

[26] J. Lohkamp, Bull. Amer. Math. Soc. (N.S.) 27 no. 2, 288 
(1992). 

[27] J. Lewandowski and T. Pawlowski, Class. Quantum 
Grav. 22, 1573 (2005). 



